**************************************
*****Code for replicating analyses in
*****Carlin, Hellwig, Love, Martinez-Gallardo, & Singer, "When does the Public Get It Right? The Information Environment and the Accuracy of Economic Sentiment" 
*****Comparative Political Studies, forthcoming
**************************************
*****code for replicating all tables and figures except Tables 2, A20, A21, & A22, and Figure 2
**************************************
*****install mipolate.pkg, grinter.pkg
**************************************
use "CHLMS_CPS_qt.dta", replace 

xtset ccode ticker
label variable yq "year"
replace opp1share=opp1seatq/totalseats
replace opp1shmc = opp1share - .3028702
gen a_unempl = q_unempl
gen a_trans = transparencyindex
gen a_fotpflip = q_fotpflip
gen a_opp1share = opp1share
recode a_unempl *=. if quarter~=4
recode a_trans *=. if quarter~=4
recode a_fotpflip *=. if quarter~=4
xtset ccode ticker 
sort ccode
by ccode: mipolate a_unempl ticker, pchip gen(c_unempl)
by ccode: mipolate a_trans ticker, pchip gen(c_trans)
by ccode: mipolate a_fotpflip ticker, pchip gen(c_fotpflip)
gen opp1pt = opp1share*100

/*CREATE STANDARDIZED MEASURES*/
gen ziccpresn = (iccpresn - 46.72)/10.15
gen zqgrowth = (qgrowth - 3.9619)/3.58357
gen zlog_q_infl = (log_q_infl - 1.44875)/.89175
gen zq_unempl = (c_unempl - 8.5621)/3.1115
gen zq_trans = (c_trans-2.002)/1.4463
gen zq_fotpflip = (c_fotpflip-59.72)/12.58163
gen zopp1share = (opp1share-.3)/.1192801
gen infoenvi = (zq_trans + zq_fotpflip + zopp1share)/3
gen objecon = (zqgrowth - zlog_q_infl - zq_unempl)/3

/*FIGURE 1*/
sort yq
label variable objecon "economic indicator index"
line iccpresn yq if year>1999 & year<2016 & ccode==160 & sample==1, lcolor(black) || line objecon yq if year>1999 & year<2016 & ccode==160 & sample==1, lcolor(gs10) yaxis(2) scheme(s1mono) title("Argentina") name(arg, replace) legend(off) nodraw
line iccpresn yq if year>1999 & year<2016 & ccode==140 & sample==1, lcolor(black) || line objecon yq if year>1999 & year<2016 & ccode==140 & sample==1, lcolor(gs10) yaxis(2) scheme(s1mono) title("Brazil") name(bra, replace) legend(off) nodraw
line iccpresn yq if year>1999 & year<2016 & ccode==155 & sample==1, lcolor(black) || line objecon yq if year>1999 & year<2016 & ccode==155 & sample==1, lcolor(gs10) yaxis(2) scheme(s1mono) title("Chile") name(chl, replace) legend(off) nodraw
line iccpresn yq if year>1999 & year<2016 & ccode==100 & sample==1, lcolor(black) || line objecon yq if year>1999 & year<2016 & ccode==100 & sample==1, lcolor(gs10) yaxis(2) scheme(s1mono) title("Colombia") name(col, replace) legend(off) nodraw
line iccpresn yq if year>1999 & year<2016 & ccode==94 & sample==1, lcolor(black) || line objecon yq if year>1999 & year<2016 & ccode==94 & sample==1, lcolor(gs10) yaxis(2) scheme(s1mono) title("Costa Rica") name(cri, replace) legend(off) nodraw
line iccpresn yq if year>1999 & year<2016 & ccode==70 & sample==1, lcolor(black) || line objecon yq if year>1999 & year<2016 & ccode==70 & sample==1, lcolor(gs10) yaxis(2) scheme(s1mono) title("Mexico") name(mex, replace) legend(off) nodraw
line iccpresn yq if year>1999 & year<2016 & ccode==135 & sample==1, lcolor(black) || line objecon yq if year>1999 & year<2016 & ccode==135 & sample==1, lcolor(gs10) yaxis(2) scheme(s1mono) title("Peru") name(per, replace) legend(off) nodraw
line iccpresn yq if year>1999 & year<2016 & ccode==92 & sample==1, lcolor(black) || line objecon yq if year>1999 & year<2016 & ccode==92 & sample==1, lcolor(gs10) yaxis(2) scheme(s1mono) title("El Salvador") name(slv, replace) legend(off) nodraw
graph combine arg bra chl col, scheme(s1mono) c(2) ycommon 
graph combine cri mex per slv, scheme(s1mono) c(2) ycommon 
	
/*TABLE 1*/
tsset ccode ticker
by ccode: gen dinfoenvi = d.infoenvi
by ccode: gen dc_trans = d.c_trans
by ccode: gen dc_fotpflip = d.c_fotpflip
by ccode: gen dopp1pt = d.opp1pt
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres infoenvi c.d.qgrowth#c.infoenvi c.l.qgrowth#c.infoenvi c.d.log_q_infl#c.infoenvi c.l.log_q_infl#c.infoenvi c.d.c_unempl#c.infoenvi c.l.c_unempl#c.infoenvi if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres c_trans c.d.qgrowth#c.c_trans c.l.qgrowth#c.c_trans c.d.log_q_infl#c.c_trans c.l.log_q_infl#c.c_trans c.d.c_unempl#c.c_trans c.l.c_unempl#c.c_trans if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres c_fotpflip c.d.qgrowth#c.c_fotpflip c.l.qgrowth#c.c_fotpflip c.d.log_q_infl#c.c_fotpflip c.l.log_q_infl#c.c_fotpflip c.d.c_unempl#c.c_fotpflip c.l.c_unempl#c.c_fotpflip if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres opp1pt c.d.qgrowth#c.opp1pt c.l.qgrowth#c.opp1pt c.d.log_q_infl#c.opp1pt c.l.log_q_infl#c.opp1pt c.d.c_unempl#c.opp1pt c.l.c_unempl#c.opp1pt if esample~=., c(psar)
	
/*FIGURES 3 & 4*/
gen diccpresn = d.iccpresn
gen dqgrowth = d.qgrowth
gen dlog_q_infl = d.log_q_infl
gen dc_unempl = d.c_unempl
gen dappqcres = d.appqcres
gen lappqcres = l.appqcres 
gen liccpresn = l.iccpresn 
gen lqgrowth = l.qgrowth 
gen llog_q_infl = l.log_q_infl
gen lc_unempl = l.c_unempl
gen dgxi = dqgrowth*infoenvi  
gen dixi = dlog_q_infl*infoenvi   
gen duxi = dc_unempl*infoenvi
gen dgxt = dqgrowth*c_trans 
gen dixt = dlog_q_infl*c_trans
gen duxt = dc_unempl*c_trans 
gen dgxf = dqgrowth*c_fotpflip  
gen dixf = dlog_q_infl*c_fotpflip  
gen duxf = dc_unempl*c_fotpflip
gen dgxo = dqgrowth*opp1pt 
gen dixo = dlog_q_infl*opp1pt  
gen duxo = dc_unempl*opp1pt
gen lgxi = lqgrowth*infoenvi  
gen lixi = llog_q_infl*infoenvi   
gen luxi = lc_unempl*infoenvi
gen lgxt = lqgrowth*c_trans 
gen lixt = llog_q_infl*c_trans
gen luxt = lc_unempl*c_trans 
gen lgxf = lqgrowth*c_fotpflip  
gen lixf = llog_q_infl*c_fotpflip  
gen luxf = lc_unempl*c_fotpflip
gen lgxo = lqgrowth*opp1pt 
gen lixo = llog_q_infl*opp1pt  
gen luxo = lc_unempl*opp1pt

/*FIGURE 3*/
label variable infoenvi "Information Environment"
quietly xtpcse diccpresn dqgrowth infoenvi ccode_dummy5 ccode_dummy17 ccode_dummy16 ccode_dummy13 ccode_dummy9 ccode_dummy7 ccode_dummy2 liccpresn lqgrowth dlog_q_infl llog_q_infl dc_unempl lc_unempl dappqcres lappqcres dgxi dixi duxi lgxi lixi luxi if esample~=., c(psar) 
grinter (dqgrowth), inter(dgxi) const(infoenvi) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Growth) graphr(c(white))
grinter (dlog_q_infl), inter(dixi) const(infoenvi) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Inflation) graphr(c(white))
grinter (dc_unempl), inter(duxi) const(infoenvi) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Unemployment) graphr(c(white))

/*FIGURE 4*/
label variable c_trans "Transparency"
quietly xtpcse diccpresn dqgrowth c_trans ccode_dummy5 ccode_dummy17 ccode_dummy16 ccode_dummy13 ccode_dummy9 ccode_dummy7 ccode_dummy2 liccpresn lqgrowth dlog_q_infl llog_q_infl dc_unempl lc_unempl dappqcres lappqcres dgxt dixt duxt lgxt lixt luxt if esample~=., c(psar) 
grinter (dqgrowth), inter(dgxt) const(c_trans) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Growth) graphr(c(white)) ylabel(-.5(.5)1.5)
grinter (dlog_q_infl), inter(dixt) const(c_trans) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Inflation) graphr(c(white)) ylabel(-9(3)9)
grinter (dc_unempl), inter(duxt) const(c_trans) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Unemployment) graphr(c(white)) ylabel(-15(5)5)
label variable c_fotpflip "Press Freedom"
quietly xtpcse diccpresn dqgrowth c_fotpflip ccode_dummy5 ccode_dummy17 ccode_dummy16 ccode_dummy13 ccode_dummy9 ccode_dummy7 ccode_dummy2 liccpresn lqgrowth dlog_q_infl llog_q_infl dc_unempl lc_unempl dappqcres lappqcres dgxf dixf duxf lgxf lixf luxf if esample~=., c(psar) 
grinter (dqgrowth), inter(dgxf) const(c_fotpflip) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Growth) graphr(c(white)) ylabel(-.5(.5)1.5)
grinter (dlog_q_infl), inter(dixf) const(c_fotpflip) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Inflation) graphr(c(white)) ylabel(-9(3)9)
grinter (dc_unempl), inter(duxf) const(c_fotpflip) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Unemployment) graphr(c(white)) ylabel(-15(5)5)
label variable opp1pt "Opposition Percent"
quietly xtpcse diccpresn dqgrowth opp1pt ccode_dummy5 ccode_dummy17 ccode_dummy16 ccode_dummy13 ccode_dummy9 ccode_dummy7 ccode_dummy2 liccpresn lqgrowth dlog_q_infl llog_q_infl dc_unempl lc_unempl dappqcres lappqcres dgxo dixo duxo lgxo lixo luxo if esample~=., c(psar) 
grinter (dqgrowth), inter(dgxo) const(opp1pt) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Growth) graphr(c(white)) ylabel(-.5(.5)1.5)
grinter (dlog_q_infl), inter(dixo) const(opp1pt) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Inflation) graphr(c(white)) ylabel(-9(3)9)
grinter (dc_unempl), inter(duxo) const(opp1pt) yline(0, lc(gs7)) lc(black) kdensity kdop(lpattern(solid) lc(gs10)) non nom nomeantext scheme(s2mono) ytitle(ME of Unemployment) graphr(c(white)) ylabel(-15(5)5)
	
	
*********************************************
****TABLES & FIGURES REPORTED IN THE SUPPLEMENTARY INFORMATION FILE
*********************************************
/*Table A1*/
sum iccpresn qgrowth log_q_infl c_unempl appqcres c_trans c_fotpflip opp1pt infoenvi if esample~=.

/*Table A2. text*/

/*Table A3*/
sort country year
merge country year using "C:\Users\thellwig\Dropbox\EAD\0 paper projects\modeling icc\CPS\replication files\Average yearly PID.dta"
sort country year quarter
drop _merge
merge country year quarter using "C:\Users\thellwig\Dropbox\EAD\0 paper projects\modeling icc\CPS\replication files\Status.dta"
drop if ccode==.
drop _merge
sort country year quarter
by country: mipolate pid_yravg year if esample~=., near gen(y_pid_rate)
order country year month quarter pid pid_yravg y_pid_rate 
xtset ccode ticker
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres infoenvi pid_yravg c.d.qgrowth#c.infoenvi c.l.qgrowth#c.infoenvi c.d.log_q_infl#c.infoenvi c.l.log_q_infl#c.infoenvi c.d.c_unempl#c.infoenvi c.l.c_unempl#c.infoenvi if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.obj l.obj d.appqcres l.appqcres infoenvi pid_yravg c.d.obj#c.infoenvi c.l.obj#c.infoenvi c.d.obj#c.pid_yravg c.l.obj#c.pid_yravg if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres infoenvi y_pid c.d.qgrowth#c.infoenvi c.l.qgrowth#c.infoenvi c.d.log_q_infl#c.infoenvi c.l.log_q_infl#c.infoenvi c.d.c_unempl#c.infoenvi c.l.c_unempl#c.infoenvi if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.obj l.obj d.appqcres l.appqcres infoenvi y_pid c.d.obj#c.infoenvi c.l.obj#c.infoenvi c.d.obj#c.y_pid c.l.obj#c.y_pid if esample~=., c(psar)

/*Table A4*/
xtunitroot ips iccpresn if esample~=., lags(aic 4)
xtunitroot ips qgrowth if esample~=., lags(aic 4)
xtunitroot ips log_q_infl if esample~=., lags(aic 4)
xtunitroot ips c_unempl if esample~=., lags(aic 4)
xtunitroot ips appqcres if esample~=., lags(aic 4)
xtunitroot ips infoenvi if esample~=., lags(aic 4)
xtunitroot ips c_trans if esample~=., lags(aic 4)
xtunitroot ips c_fotpflip if esample~=., lags(aic 4)
xtunitroot ips opp1pt if esample~=., lags(aic 4)

/*Table A5*/
gen lincomb1 = -.216*l.iccpresn + .372*d.qgrowth + .281*l.qgrowth - 2.457*d.log_q_infl - .634*l.log_q_infl -3.12*d.c_unempl -.232*l.c_unempl +.089*d.appqcres +.033*l.appqcres 
gen lincomb2 = -.272*l.iccpresn + .284*d.qgrowth + .326*l.qgrowth + .902*d.log_q_infl + .601*l.log_q_infl -.314*d.c_unempl -.517*l.c_unempl +.098*d.appqcres +.067*l.appqcres + 4.174*infoenvi - .445*c.d.qgrowth#c.infoenvi + .051*c.l.qgrowth#c.infoenvi - 4.945*c.d.log_q_infl#c.infoenvi -1.825*c.l.log_q_infl#c.infoenvi -8.989*c.d.c_unempl#c.infoenvi -.027*c.l.c_unempl#c.infoenvi
gen lincomb3 = -.194*l.iccpresn + .116*d.qgrowth + .048*l.qgrowth + 4.033*d.log_q_infl + 1.124*l.log_q_infl -4.42*d.c_unempl +.32*l.c_unempl +.087*d.appqcres +.034*l.appqcres + 1.34*c_trans + .061*c.d.qgrowth#c.c_trans + .061*c.l.qgrowth#c.c_trans - 2.006*c.d.log_q_infl#c.c_trans -.508*c.l.log_q_infl#c.c_trans + .376*c.d.c_unempl#c.c_trans -.1*c.l.c_unempl#c.c_trans
gen lincomb4 = -.226*l.iccpresn -1.03*d.qgrowth -.254*l.qgrowth + 11.775*d.log_q_infl + 3.155*l.log_q_infl -.096*d.c_unempl +.309*l.c_unempl +.089*d.appqcres +.041*l.appqcres + .12*c_fotpflip + .026*c.d.qgrowth#c.c_fotpflip + .01*c.l.qgrowth#c.c_fotpflip - .219*c.d.log_q_infl#c.c_fotpflip -.056*c.l.log_q_infl#c.c_fotpflip -.046*c.d.c_unempl#c.c_fotpflip +.0003*c.l.c_unempl#c.c_fotpflip
gen lincomb5 = -.242*l.iccpresn + .127*d.qgrowth + .039*l.qgrowth + 2.648*d.log_q_infl + .268*l.log_q_infl + 3.347*d.c_unempl - .28*l.c_unempl +.096*d.appqcres +.044*l.appqcres + .045*opp1pt + .007*c.d.qgrowth#c.opp1pt +.007*c.l.qgrowth#c.opp1pt - .137*c.d.log_q_infl#c.opp1pt -.026*c.l.log_q_infl#c.opp1pt -.311*c.d.c_unempl#c.opp1pt -.00003*c.l.c_unempl#c.opp1pt
xtunitroot ips lincomb1 if esample~=., lags(aic 4)
xtunitroot ips lincomb2 if esample~=., lags(aic 4)
xtunitroot ips lincomb3 if esample~=., lags(aic 4)
xtunitroot ips lincomb4 if esample~=., lags(aic 4)
xtunitroot ips lincomb5 if esample~=., lags(aic 4)

/*Table A6 - checking residuals*/
qui xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres if esample~=., c(psar)
predict fittedM1
gen residM1 = d.iccpresn - fittedM1
label variable residM1 "T1 M1 residuals"
qui xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres infoenvi c.d.qgrowth#c.infoenvi c.l.qgrowth#c.infoenvi c.d.log_q_infl#c.infoenvi c.l.log_q_infl#c.infoenvi c.d.c_unempl#c.infoenvi c.l.c_unempl#c.infoenvi if esample~=., c(psar)
predict fittedM2
gen residM2 = d.iccpresn - fittedM2
label variable residM2 "T1 M2 residuals"
qui xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres c_trans c.d.qgrowth#c.c_trans c.l.qgrowth#c.c_trans c.d.log_q_infl#c.c_trans c.l.log_q_infl#c.c_trans c.d.c_unempl#c.c_trans c.l.c_unempl#c.c_trans if esample~=., c(psar)
predict fittedM3
gen residM3 = d.iccpresn - fittedM3
label variable residM3 "T1 M3 residuals"
qui xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres c_fotpflip c.d.qgrowth#c.c_fotpflip c.l.qgrowth#c.c_fotpflip c.d.log_q_infl#c.c_fotpflip c.l.log_q_infl#c.c_fotpflip c.d.c_unempl#c.c_fotpflip c.l.c_unempl#c.c_fotpflip if esample~=., c(psar)
predict fittedM4
gen residM4 = d.iccpresn - fittedM4
label variable residM4 "T1 M4 residuals"
qui xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres opp1pt c.d.qgrowth#c.opp1pt c.l.qgrowth#c.opp1pt c.d.log_q_infl#c.opp1pt c.l.log_q_infl#c.opp1pt c.d.c_unempl#c.opp1pt c.l.c_unempl#c.opp1pt if esample~=., c(psar)
predict fittedM5
gen residM5 = d.iccpresn - fittedM5
label variable residM5 "T1 M5 residuals"
xtunitroot ips residM1 if esample~=., lags(aic 4)
xtunitroot ips residM2 if esample~=., lags(aic 4)
xtunitroot ips residM3 if esample~=., lags(aic 4)
xtunitroot ips residM4 if esample~=., lags(aic 4)
xtunitroot ips residM5 if esample~=., lags(aic 4)

/*Table A7 - Residuals one country at a time*/
dfuller residM1 if ccode==160 & esample~=., lags(1)
dfuller residM1 if ccode==140 & esample~=., lags(1)
dfuller residM1 if ccode==155 & esample~=., lags(1)
dfuller residM1 if ccode==100 & esample~=., lags(1)
dfuller residM1 if ccode==94 & esample~=., lags(1)
dfuller residM1 if ccode==70 & esample~=., lags(1)
dfuller residM1 if ccode==135 & esample~=., lags(1)
dfuller residM1 if ccode==92 & esample~=., lags(1)

/*Figure A1*/
line residM1 yq if ccode==160 & esample~=., scheme(s1mono) title("Argentina") name(arg_res1, replace) legend(off) nodraw
line residM1 yq if ccode==140 & esample~=., scheme(s1mono) title("Brazil") name(bra_res1, replace) legend(off) nodraw
line residM1 yq if ccode==155 & esample~=., scheme(s1mono) title("Chile") name(chl_res1, replace) legend(off) nodraw
line residM1 yq if ccode==100 & esample~=., scheme(s1mono) title("Colombia") name(col_res1, replace) legend(off) nodraw
line residM1 yq if ccode==94 & esample~=., scheme(s1mono) title("Costa Rica") name(cri_res1, replace) legend(off) nodraw
line residM1 yq if ccode==70 & esample~=., scheme(s1mono) title("Mexico") name(mex_res1, replace) legend(off) nodraw
line residM1 yq if ccode==135 & esample~=., scheme(s1mono) title("Peru") name(per_res1, replace) legend(off) nodraw
line residM1 yq if ccode==92 & esample~=., scheme(s1mono) title("El Salvador") name(slv_res1, replace) legend(off) nodraw
graph combine arg_res1 bra_res1 chl_res1 col_res1 cri_res1 mex_res1 per_res1 slv_res1, scheme(s1mono) ycommon

/*Table A8*/
sum dinfoenvi dc_trans dc_fotpflip dopp1pt if esample~=.

/*Table A9*/
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.infoenvi l.infoenvi c.d.qgrowth#c.l.infoenvi c.l.qgrowth#c.d.infoenvi c.d.qgrowth#c.d.infoenvi c.l.qgrowth#c.l.infoenvi if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.infoenvi l.infoenvi c.d.log_q_infl#c.l.infoenvi c.l.log_q_infl#c.d.infoenvi c.d.log_q_infl#c.d.infoenvi c.l.log_q_infl#c.l.infoenvi if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.infoenvi l.infoenvi c.d.c_unempl#c.l.infoenvi c.l.c_unempl#c.d.infoenvi c.d.c_unempl#c.d.infoenvi c.l.c_unempl#c.l.infoenvi if esample~=., c(psar)

/*Table A10*/	
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.c_trans l.c_trans c.d.qgrowth#c.l.c_trans c.l.qgrowth#c.d.c_trans c.d.qgrowth#c.d.c_trans c.l.qgrowth#c.l.c_trans if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.c_trans l.c_trans c.d.log_q_infl#c.l.c_trans c.l.log_q_infl#c.d.c_trans c.d.log_q_infl#c.d.c_trans c.l.log_q_infl#c.l.c_trans if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.c_trans  l.c_trans c.d.c_unempl#c.l.c_trans c.l.c_unempl#c.d.c_trans c.d.c_unempl#c.d.c_trans c.l.c_unempl#c.l.c_trans if esample~=., c(psar)

/*Table A11*/	
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.c_fotpf l.c_fotpf c.d.qgrowth#c.l.c_fotpf c.l.qgrowth#c.d.c_fotpf c.d.qgrowth#c.d.c_fotpf c.l.qgrowth#c.l.c_fotpf if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.c_fotpf  l.c_fotpf c.d.log_q_infl#c.l.c_fotpf c.l.log_q_infl#c.d.c_fotpf c.d.log_q_infl#c.d.c_fotpf c.l.log_q_infl#c.l.c_fotpf if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.c_fotpf  l.c_fotpf  c.d.c_unempl#c.l.c_fotpf c.l.c_unempl#c.d.c_fotpf c.d.c_unempl#c.d.c_fotpf c.l.c_unempl#c.l.c_fotpf if esample~=., c(psar)

/*Table A12*/	
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.opp1pt l.opp1pt c.d.qgrowth#c.l.opp1pt c.l.qgrowth#c.d.opp1pt c.d.qgrowth#c.d.opp1pt c.l.qgrowth#c.l.opp1pt if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.opp1pt l.opp1pt c.d.log_q_infl#c.l.opp1pt c.l.log_q_infl#c.d.opp1pt c.d.log_q_infl#c.d.opp1pt c.l.log_q_infl#c.l.opp1pt if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres d.opp1pt l.opp1pt c.d.c_unempl#c.l.opp1pt c.l.c_unempl#c.d.opp1pt c.d.c_unempl#c.d.opp1pt c.l.c_unempl#c.l.opp1pt if esample~=., c(psar)

/**Table A13*/
save "CHLMS_CPS_qtm.dta", replace 
/*collapse data into means*/
collapse iccpresn qgrowth log_q_infl c_unempl appqcres infoenvi c_trans c_fotpflip opp1pt esample, by(ccode year)
save "CHLMS_CPS_qtms.dta", replace 
tsset ccode year
reg d.iccpresn i.ccode d.qgrowth d.log_q_infl d.c_unempl d.appqcres if esample~=.
reg d.iccpresn i.ccode d.qgrowth d.log_q_infl d.c_unempl d.appqcres infoenvi c.d.qgrowth#c.infoenvi c.d.log_q_infl#c.infoenvi c.d.c_unempl#c.infoenvi if esample~=.
reg d.iccpresn i.ccode d.qgrowth d.log_q_infl d.c_unempl d.appqcres c_trans c.d.qgrowth#c.c_trans c.d.log_q_infl#c.c_trans c.d.c_unempl#c.c_trans if esample~=.
reg d.iccpresn i.ccode d.qgrowth d.log_q_infl d.c_unempl d.appqcres c_fotpflip c.d.qgrowth#c.c_fotpflip c.d.log_q_infl#c.c_fotpflip c.d.c_unempl#c.c_fotpflip if esample~=.
reg d.iccpresn i.ccode d.qgrowth d.log_q_infl d.c_unempl d.appqcres opp1pt c.d.qgrowth#c.opp1pt c.d.log_q_infl#c.opp1pt c.d.c_unempl#c.opp1pt if esample~=.
save "CHLMS_CPS_qtms.dta", replace 
clear
use "CHLMS_CPS_qtm.dta", replace 

/*Table A14*/
xtpcse d.iccpresn l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres if esample~=., c(psar)
xtpcse d.iccpresn l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres infoenvi c.d.qgrowth#c.infoenvi c.l.qgrowth#c.infoenvi c.d.log_q_infl#c.infoenvi c.l.log_q_infl#c.infoenvi c.d.c_unempl#c.infoenvi c.l.c_unempl#c.infoenvi if esample~=., c(psar)
xtpcse d.iccpresn l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres c_trans c.d.qgrowth#c.c_trans c.l.qgrowth#c.c_trans c.d.log_q_infl#c.c_trans c.l.log_q_infl#c.c_trans c.d.c_unempl#c.c_trans c.l.c_unempl#c.c_trans if esample~=., c(psar)
xtpcse d.iccpresn l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres c_fotpflip c.d.qgrowth#c.c_fotpflip c.l.qgrowth#c.c_fotpflip c.d.log_q_infl#c.c_fotpflip c.l.log_q_infl#c.c_fotpflip c.d.c_unempl#c.c_fotpflip c.l.c_unempl#c.c_fotpflip if esample~=., c(psar)
xtpcse d.iccpresn l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres opp1pt c.d.qgrowth#c.opp1pt c.l.qgrowth#c.opp1pt c.d.log_q_infl#c.opp1pt c.l.log_q_infl#c.opp1pt c.d.c_unempl#c.opp1pt c.l.c_unempl#c.opp1pt if esample~=., c(psar)

/*Table A15*/
gen q_fotpaflip = 100 - q_fotpa_adj
gen q_fotpbflip = 100 - q_fotpb_adj
gen q_fotpcflip = 100 - q_fotpc_adj
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres q_fotpaflip c.d.qgrowth#c.q_fotpaflip c.l.qgrowth#c.q_fotpaflip c.d.log_q_infl#c.q_fotpaflip c.l.log_q_infl#c.q_fotpaflip c.d.c_unempl#c.q_fotpaflip c.l.c_unempl#c.q_fotpaflip if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres q_fotpbflip c.d.qgrowth#c.q_fotpbflip c.l.qgrowth#c.q_fotpbflip c.d.log_q_infl#c.q_fotpbflip c.l.log_q_infl#c.q_fotpbflip c.d.c_unempl#c.q_fotpbflip c.l.c_unempl#c.q_fotpbflip if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres q_fotpcflip c.d.qgrowth#c.q_fotpcflip c.l.qgrowth#c.q_fotpcflip c.d.log_q_infl#c.q_fotpcflip c.l.log_q_infl#c.q_fotpcflip c.d.c_unempl#c.q_fotpcflip c.l.c_unempl#c.q_fotpcflip if esample~=., c(psar)

/*Table A16*/ 
gen q_w_itimc = q_w_iti - 70.2
gen q_w_quantmc = q_w_quant - 77.4
gen q_w_qualcmc = q_w_qualc - 68.4
pwcorr infoenvi q_w_iti q_w_quant q_w_qualc if esample~=., sig obs

/*Table A17*/
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres q_w_itimc c.d.qgrowth#c.q_w_itimc c.l.qgrowth#c.q_w_itimc c.d.log_q_infl#c.q_w_itimc c.l.log_q_infl#c.q_w_itimc c.d.c_unempl#c.q_w_itimc c.l.c_unempl#c.q_w_itimc if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres q_w_quantmc c.d.qgrowth#c.q_w_quantmc c.l.qgrowth#c.q_w_quantmc c.d.log_q_infl#c.q_w_quantmc c.l.log_q_infl#c.q_w_quantmc c.d.c_unempl#c.q_w_quantmc c.l.c_unempl#c.q_w_quantmc if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres q_w_qualcmc c.d.qgrowth#c.q_w_qualcmc c.l.qgrowth#c.q_w_qualcmc c.d.log_q_infl#c.q_w_qualcmc c.l.log_q_infl#c.q_w_qualcmc c.d.c_unempl#c.q_w_qualcmc c.l.c_unempl#c.q_w_qualcmc if esample~=., c(psar)
/*extracting qual from infoenvi*/
reg infoenvi q_w_qualcmc
predict ie_noqualc if e(sample), res
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres ie_noqualc c.d.qgrowth#c.ie_noqualc c.l.qgrowth#c.ie_noqualc c.d.log_q_infl#c.ie_noqualc c.l.log_q_infl#c.ie_noqualc c.d.c_unempl#c.ie_noqualc c.l.c_unempl#c.ie_noqualc if esample~=., c(psar)

/*Table A18*/
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres q_statcap infoenvi c.d.qgrowth#c.infoenvi c.d.log_q_infl#c.infoenvi c.d.c_unempl#c.infoenvi if esample~=., c(psar)
predict sqsample2 if e(sample)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres q_statcap infoenvi c.d.qgrowth#c.q_statcap c.d.log_q_infl#c.q_statcap c.d.c_unempl#c.q_statcap if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres q_statcap infoenvi c.d.qgrowth#c.infoenvi c.d.log_q_infl#c.infoenvi c.d.c_unempl#c.infoenvi c.d.qgrowth#c.q_statcap c.d.log_q_infl#c.q_statcap c.d.c_unempl#c.q_statcap if esample~=., c(psar)

/*Table A19*/
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres infoenvi c.d.qgrowth#c.infoenvi c.l.qgrowth#c.infoenvi c.d.log_q_infl#c.infoenvi c.l.log_q_infl#c.infoenvi c.d.c_unempl#c.infoenvi c.l.c_unempl#c.infoenvi govstab c.d.qgrowth#c.govstab c.l.qgrowth#c.govstab c.d.log_q_infl#c.govstab c.l.log_q_infl#c.govstab c.d.c_unempl#c.govstab c.l.c_unempl#c.govstab if esample~=., c(psar)
xtpcse d.iccpresn i.ccode l.iccpresn d.qgrowth l.qgrowth d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres infoenvi c.d.qgrowth#c.infoenvi c.l.qgrowth#c.infoenvi c.d.log_q_infl#c.infoenvi c.l.log_q_infl#c.infoenvi c.d.c_unempl#c.infoenvi c.l.c_unempl#c.infoenvi ers c.d.qgrowth#c.ers c.l.qgrowth#c.ers c.d.log_q_infl#c.ers c.l.log_q_infl#c.ers c.d.c_unempl#c.ers c.l.c_unempl#c.ers if esample~=., c(psar)


/**********Figure A2*************/
/*LRE INFOENVI GROWTH*/
xtpcse d.iccpresn d.qgrowth l.qgrowth infoenvi c.d.qgrowth#c.infoenvi c.l.qgrowth#c.infoenvi l.iccpresn d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres c.d.log_q_infl#c.infoenvi c.l.log_q_infl#c.infoenvi c.d.c_unempl#c.infoenvi c.l.c_unempl#c.infoenvi i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist -.75(.25)1.25 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Information Environment") ytitle("LRE of Growth") legend(off) scheme(s1mono) name(lre_grie, replace) nodraw

/*LRE INFOENVI INFLATION*/
xtpcse d.iccpresn d.log_q_infl l.log_q_infl infoenvi c.d.log_q_infl#c.infoenvi c.l.log_q_infl#c.infoenvi l.iccpresn c.d.qgrowth#c.infoenvi c.l.qgrowth#c.infoenvi d.qgrowth l.qgrowth d.c_unempl l.c_unempl d.appqcres l.appqcres c.d.c_unempl#c.infoenvi c.l.c_unempl#c.infoenvi i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist -.75(.25)1.25 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Information Environment") ytitle("LRE of Inflation") legend(off) scheme(s1mono) name(lre_inie, replace) nodraw

/*LRE INFOENVI UNEMPLOYMENT*/
xtpcse d.iccpresn d.c_unempl l.c_unempl infoenvi c.d.c_unempl#c.infoenvi c.l.c_unempl#c.infoenvi l.iccpresn d.qgrowth l.qgrowth c.d.qgrowth#c.infoenvi c.l.qgrowth#c.infoenvi d.log_q_infl l.log_q_infl d.appqcres l.appqcres c.d.log_q_infl#c.infoenvi c.l.log_q_infl#c.infoenvi  i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist -.75(.25)1.25 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Information Environment") ytitle("LRE of Unemployment") legend(off) scheme(s1mono) name(lre_unie, replace) nodraw

/*LRE TRANSPARENCY GROWTH*/
xtpcse d.iccpresn d.qgrowth l.qgrowth c_trans c.d.qgrowth#c.c_trans c.l.qgrowth#c.c_trans l.iccpresn d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres c.d.log_q_infl#c.c_trans c.l.log_q_infl#c.c_trans c.d.c_unempl#c.c_trans c.l.c_unempl#c.c_trans i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist 1.5(.5)5.5 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Transparency") ytitle("LRE of Growth") legend(off) scheme(s1mono) name(lre_grtr, replace) nodraw

/*LRE TRANSPARENCY INFLATION*/
xtpcse d.iccpresn d.log_q_infl l.log_q_infl c_trans c.d.log_q_infl#c.c_trans c.l.log_q_infl#c.c_trans l.iccpresn d.qgrowth l.qgrowth d.c_unempl l.c_unempl c.d.qgrowth#c.c_trans c.l.qgrowth#c.c_trans  d.appqcres l.appqcres c.d.c_unempl#c.c_trans c.l.c_unempl#c.c_trans i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist 1.5(.5)5.5 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Transparency") ytitle("LRE of Inflation") legend(off) scheme(s1mono) name(lre_intr, replace) nodraw

/*LRE TRANSPARENCY UNEMPLOYMENT*/
xtpcse d.iccpresn d.c_unempl l.c_unempl c_trans c.d.log_q_infl#c.c_trans c.l.log_q_infl#c.c_trans l.iccpresn d.log_q_infl l.log_q_infl d.appqcres l.appqcres d.qgrowth l.qgrowth c.d.qgrowth#c.c_trans c.l.qgrowth#c.c_trans c.d.c_unempl#c.c_trans c.l.c_unempl#c.c_trans i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist 1.5(.5)5.5 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Transparency") ytitle("LRE of Unemployment") legend(off) scheme(s1mono) name(lre_untr, replace) nodraw

/*LRE PRESS GROWTH*/
xtpcse d.iccpresn d.qgrowth l.qgrowth c_fotpflip c.d.qgrowth#c.c_fotpflip c.l.qgrowth#c.c_fotpflip l.iccpresn d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres c.d.log_q_infl#c.c_fotpflip c.l.log_q_infl#c.c_fotpflip c.d.c_unempl#c.c_fotpflip c.l.c_unempl#c.c_fotpflip i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist 40(5)80 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Press Freedom") ytitle("LRE of Growth") legend(off) scheme(s1mono) name(lre_grfp, replace) nodraw

/*LRE PRESS INFLATION*/
xtpcse d.iccpresn d.log_q_infl l.log_q_infl c_fotpflip c.d.log_q_infl#c.c_fotpflip c.l.log_q_infl#c.c_fotpflip l.iccpresn d.qgrowth l.qgrowth c.d.qgrowth#c.c_fotpflip c.l.qgrowth#c.c_fotpflip d.c_unempl l.c_unempl d.appqcres l.appqcres c.d.c_unempl#c.c_fotpflip c.l.c_unempl#c.c_fotpflip i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist 40(5)80 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Press Freedom") ytitle("LRE of Inflation") legend(off) scheme(s1mono) name(lre_infp, replace) nodraw

/*LRE PRESS UNEMPLOYMENT*/
xtpcse d.iccpresn d.c_unempl l.c_unempl c_fotpflip c.d.c_unempl#c.c_fotpflip c.l.c_unempl#c.c_fotpflip l.iccpresn d.qgrowth l.qgrowth c.d.qgrowth#c.c_fotpflip c.l.qgrowth#c.c_fotpflip d.appqcres l.appqcres  d.log_q_infl l.log_q_infl c.d.log_q_infl#c.c_fotpflip c.l.log_q_infl#c.c_fotpflip i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist 40(5)80 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Press Freedom") ytitle("LRE of Unemployment") legend(off) scheme(s1mono) name(lre_unfp, replace) nodraw

/*LRE OPPOSITION GROWTH*/
xtpcse d.iccpresn d.qgrowth l.qgrowth opp1pt c.d.qgrowth#c.opp1pt c.l.qgrowth#c.opp1pt l.iccpresn d.log_q_infl l.log_q_infl d.c_unempl l.c_unempl d.appqcres l.appqcres c.d.log_q_infl#c.opp1pt c.l.log_q_infl#c.opp1pt c.d.c_unempl#c.opp1pt c.l.c_unempl#c.opp1pt i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist 5(5)45 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Opposition Percent") ytitle("LRE of Growth") legend(off) scheme(s1mono) name(lre_grop, replace) nodraw

/*LRE OPPOSITION PERCENT INFLATION*/
xtpcse d.iccpresn d.log_q_infl l.log_q_infl opp1pt c.d.log_q_infl#c.opp1pt c.l.log_q_infl#c.opp1pt l.iccpresn d.qgrowth l.qgrowth c.d.qgrowth#c.opp1pt c.l.qgrowth#c.opp1pt d.c_unempl l.c_unempl d.appqcres l.appqcres c.d.c_unempl#c.opp1pt c.l.c_unempl#c.opp1pt i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist 5(5)45 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Opposition Percent") ytitle("LRE of Inflation") legend(off) scheme(s1mono) name(lre_inop, replace) nodraw

/*LRE OPPOSITION UNEMPLOYMENT*/
xtpcse d.iccpresn d.c_unempl l.c_unempl opp1pt c.d.c_unempl#c.opp1pt c.l.c_unempl#c.opp1pt  l.iccpresn d.qgrowth l.qgrowth c.d.qgrowth#c.opp1pt c.l.qgrowth#c.opp1pt d.log_q_infl l.log_q_infl d.appqcres l.appqcres c.d.log_q_infl#c.opp1pt c.l.log_q_infl#c.opp1pt i.ccode if esample~=., c(psar)

matrix b=e(b) 
matrix V=e(V)
scalar b1=b[1,1]
scalar b2=b[1,2]
scalar b3=b[1,3]
scalar b4=b[1,4]
scalar b5=b[1,5]
scalar b6=b[1,6]
scalar varb1=V[1,1]
scalar varb2=V[2,2]
scalar varb4=V[4,4]
scalar varb5=V[5,5]
scalar varb6=V[6,6]
scalar covb1b4=V[1,4] 
scalar covb2b5=V[2,5]

capture drop vectaxis me lower_ci upper_ci
gen vectaxis = .
gen me = .
gen lower_ci = .
gen upper_ci = .
local b = 1

foreach i of numlist 5(5)45 {
	local n = 1
	qui gen conb = (b2+(b5*(`i')))/(-1*b6)
	local conse = sqrt((1/(b6^2))*(varb2+varb5*((`i')^2)+2*covb2b5*(`i')) + (((b2+(b5*(`i')))^2)/b6^4)*varb6 - 2*(b2+(b5*(`i')))/(-1*b6^3)*(varb2+varb5*((`i')^2)+2*covb2b5*(`i'))*varb6)
	qui gen x=1.96*`conse'
 	qui gen upper=conb+x
	qui gen lower=conb-x
	qui sum conb
		local me = `r(mean)'
		replace me = `me' in `b'
	qui sum lower 
		local lower_ci = `r(mean)'
		replace lower_ci = `lower_ci' in `b'
	qui sum upper
		local upper_ci = `r(mean)'
		replace upper_ci = `upper_ci' in `b'
	di
	di
	di "The marginal effect for consumer confidence is " `me'
	di "The confidence interval is " `lower_ci' 	"	and		"	`upper_ci'
	di "The modifying variable is = 	" `i'
	di 
	di
	capture drop conb* x* upper lower
	replace vectaxis = `i' in `b'
	local b = `b' + 1
	local n = `n' + 1
}

twoway (line me vectaxis, clpattern(solid) sort) (line lower_ci vectaxis, clpattern(dash) sort) (line upper_ci vectaxis, clpattern(dash) sort), yline(0) xtitle("Opposition Percent") ytitle("LRE of Unemployment") legend(off) scheme(s1mono) name(lre_unop, replace) nodraw

graph combine lre_grie lre_inie lre_unie lre_grtr lre_intr lre_untr lre_grfp lre_infp lre_unfp lre_grop lre_inop lre_unop, c(3) scheme(s1mono)
